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Abstract 



bJ)| We demonstrate the effects of an induced disorder (or a free-orientation : 6 which is 

^ ■ related to the relative direction of scattering of particles w.r.t. to the normal of the propa- 

\ gating plane-wave front) upon the possible resonance of the plane (sound) wave propagating 

in Bose gases by using the quantum kinetic equations. We firstly present the diverse dis- 
, persion relations obtained by the relevant Pauli-blocking parameter B (which describes the 

Bose particles when B is positive) and the free-orientation 6 and then, based on the acoustic 
^ , analog, address the possible resonant states. 

! PACS : 42.25.Bs; 02.30.Jr; 05.30.Jp; 34.50.-s; 42.50.Ar; 42.50.Gy; 43.35.Ae; 43.35. Gk; 

O ■ 67.40.Mj 

oo 
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O ■ 1 Introduction 



Oh! Bose-Einstein condensation, as is known, occurs when a macroscopic number of bosons piles 

down to the same lowest single-particle state. In the presence of a trapping potential, there 
arises a whole spectrum of energy levels {Ei) associated with the stationary solutions to the 



Gross-Pitaevskii equation. Each of these stationary solutions is, by definition, a topological 
coherent mode [1-2]. The lowest energy level corresponds to the standard Bose-Einstein conden- 
^ ' sate, while the higher states describe various nonground-state condensates. The latter can be 

' generated by alternating fields, whose frequencies should be in resonance with the related tran- 

sition frequencies [2-3] . It is feasible to realize an oscillatory modulation of either the trapping 
potential or of the atomic scattering length. The trapping potential can be either single-well or 
multi-well. Meanwhile, the theoretically predicted low-energy collective oscillations of the Bose 
condensate have been experimentally confirmed by laser imaging techniques [1]. Moreover, the 
dynamics of the collective oscillations of the condensate has been theoretically studied beyond 
the linear regime, showing strong enhancement of the amplitude dependence of frequencies in 
presence of resonances (cf. Pitaevskii and Stringari in [1]). 

Very recently, optical lattices and atomic transport therein have attracted new attention with 
the achievement of Bose-Einstein condensation by purely optical means and the observation of 
a superfiuid to Mott insulator transition in a gas of ultracold atoms [4]. However, no detailed 
study of the dependence of spatial diffusion on the different directions of an anisotropic lattice 
have been performed so far. In fact, cold atomic samples also constitute an ideal system for 
the study of complex nonlinear phenomena and allows for the direct observation of wave packet 
dynamics in real space on a macroscopic scale [1-4]. Above mentioned progresses are closely 
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relevant to the quantum-mechanic many-body phenomena subjected to extremely confinements. 
One special interest is related to the dynamic resonances in cold Bose gases [2-3] . 
On the other hand, emerging interests in the wave propagation in the random, disordered, and 
granular media under the influence of spatial confinement as well as studies of collision phe- 
nomena in rather cold gases, e.g., Bose gases have stimulated intensive researches recently [5-8]. 
For instance, Varshni investigated the spectra of helium (bosonic particles) at high pressures 
[5] . The properties of helium atoms confined to move in restricted geometry have been of con- 
siderable interest during the last several decades [9-10]. The confined helium atom, however, 
needs much more complicated calculations than the confined hydrogen atom. Varshni consid- 
ered the effect of confinement on some of the lower lying excited states and the resulting optical 
spectrum from transitions between these states. His results could be applicable to (a) bubbles 
of helium implanted in a variety of materials (say, metals [11]), the physical state of helium in 
these bubbles had to be inferred by indirect means; (b) high-pressure helium plasmas [12] (to 
provide a quantitative diagnostic for plasma density); (c) astrophysics [5]. 

We noticed that acoustical analogs [13-14] considering continuum-mechanic and quantum-mechanic 
approaches are currently in rapid progress for both theories and measurements. The energy (E) 
adopted in quantum-mechanic formulation directly links to the acoustical frequency (w) consid- 
ered in classic (continuum-mechanic) or semi-classic (kinetic) formulation due to the existence 
of an acoustical analogy [14] . Classical systems could be used to study time-dependent potential 
fields and nonlinear effects, which are very difficult and time-consuming to treat numerically or 
analytically in quantum-mechanic ways. Motivated by the need to understand the wave dynam- 
ics of cold Bose gases [15-16] at the finite temperature (T > Tc) under strong confinement (which 
is beyond the hard-core approximation we adopted here although our approach can provide the 
theoretical understanding about how the collision frequency or the rarefaction parameter can 
tune the wave dispersion), an investigation for observing the microscopic localization (which 
will induce resonances in confined microdomain) [17-20] using the quantum kinetic model was 
performed and will be presented here by taking into account the acoustic analog [14]. 
In this paper, considering the quantum (discrete) kinetic model and the Uehling-Uhlenbeck col- 
lision term which could describe the collision of a gas of hard-sphere Bose-particles by tuning a 
Pauil-bloeking parameter 7 [21] (via a Pauli blocking factor of the form 1 -|- 7A''o with A'^o being 
a normalized number density giving the number of particles per cell in phase space), we will 
investigate the possible resonant states when plane waves propagates in (hard-sphere) Bose gases 
by introducing a disorder (say, induced by a high pressure or external field) or free-orientation 
{0 which is related to the relative direction of scattering of particles with respect to the normal 
of the propagating plane-wave front) and then obtaining the diverse dispersion relations which 
can thus be applied to the acoustical analogs [13-14]). This presentation will provide more clues 
to the studies of the quantum wave dynamics in Bose gases under suitable confined conditions 
and the possible appearance of the resonant states which are linked to the particles (number) 
density, induced disorder or relative free-orientation {0) and their energy states (the analogy 
between E and uj) [2-3] . 

The necessary verification of our approach with the previous available approaches [22-27] will 
be checked firstly. Our preliminary results show that, for the dispersion part (ratio of the 
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(phase) wave speed to that of hydrodynamical hmit), the quahtative agreement (for hydrody- 
namic regime) with Cowell et aZ.'s result [26] or Andrews et al.'s result [22] is quite good. As for 
the damping or attenuation part of ours, the qualitative agreement with Jackson and Zaremba's 
(m=0 mode) result [28] looks also rather good. Our results show that, as illustrated below, the 
localized or resonant states will strongly depend on the energy {E or uj), the effective scattering 
cross-section, the number density once the disorder and Pauli-blocking parameters are selected. 
We also found that, as 6 = tt/A, 9 being a disorder parameter, there exist possible resonant 
states which are similar to those reported in [13-14,17-19,29]. 

2 Formulations 

We firstly make the following assumptions before we introduce the general equations of our 
model [14,25] : 

(1) Consider a gas of identical particles of unit mass and a shape of a disk of diameter d, then 
each particle i, i = 1, - ■ ■ ,N, is characterized by the position of its center qi and its velocity Vi. 
We also have the geometric limitations : \qi — qj\ > d, i ^ j . This is illustrated schematically in 
Fig. 1. 

(2) Each particle moves in the plane with a relative velocity (in a centre-of-mass coordinate 
system) belonging to a discrete set V of 4 velocities with only one speed (due to conservation of 
momentum and energy) in the plane (4 possible different directions) during a binary encounter. 
The velocity modulus c is a reference speed depending on the reference frame and specific dis- 
tribution of particles. 

(3) The coUisional mechanism is that of rigid spheres, that is, the particles scatter elastically and 
they change their phase states instantaneously, preserving momentum. Only binary collisions 
are considered, since multiple collisions here arc negligible. 

The collisions between two particles (say i and j) take place when they are located at Qi and 
Qj = Qi — dn, where n is the unit vector joining their centers. After collisions the particles scatter, 
preserving momentum, in the directions allowed by the discrete set V. In other words, particles 
change according to {qi,Vi) {qi,v*), (qj,Vj) {Qj,Vj) . The collision is uniquely determined 
if the incoming velocity and the impact angle V') V' ^ [~^/2,7i"/2], are known, which is defined as 
the angle between Vi and n or n(^) = (cos [tp + {k — l)7r/2], sin + {k — l)7r/2]). A; = 1, • • • , 4. 
Prom the selected velocities we have two classes of encounters, i.e. {vi,Vj) = and {vi,Vj) = 
—cP, respectively. 

(a) In the first class momentum conservation implies only : encounters at 7r/2 with exchange of 
velocities Vi = ^ v* = v'''^^, Vj = v'''^^ Vj = v'' , k = 1, ■ ■ ■ , 4, m the case ip € [— 7r/2, 0], 
and Vi = v^ ^ v* = v^'^^, vj = ■u'^+^ Vj =v^,m the case V' € [0, 7r/2]. 

(b) Similarly, {vi,Vj) = —c?\ 

(i) Head-on encounters with impact angle ip = {) such that Vi = ^ v* = i;*^"'"^, Vj = v'^'^^ 
V* = v'', k = 1, - ■ ■ ,4, 

(ii) Head-on cncoTinters with impact angle tJj ^ such that vi = ^ v* = v^^'^, Vj = 

V* = v^+^, if ^ e [-vr/2, ^vi = v^ ^ v* = v''+'-^, Vj = v^^"^ v* = v''+'^, if ^ € [0, 7r/2]. 
The schematic presentation is shown in Fig 2. For grazing collisions, that is {n,Vi)= {n,Vj) = 
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0, we put V* = Vi, V* = Vj. 




Fig. 1 Schematic diagram of a collision Fig. 2 A head-on collision 

We then assume that the gas (i.e., only a two-body encounter is possible) is composed of identical 
hard-sphere particles of the same mass [14,25]. The possible velocities of these (N) particles are 
restricted to, e.g., : Ui, U2, • • • , Up, p is a finite positive integer (p 7^ N). That is to say, only 
the velocity space is discretized, the space and time variables are still continuous. The discrete 
number density of particles are denoted by 7Vj(x, associated with the velocity Uj at point x 
and time t. If only nonlinear binary collisions and the evolution of Ni are considered, we have 

^+Ui-VNi = Fi^j2T.(^kiNkNi-A^^NiNj), i = l, (1) 

where and {k,l) arc admissible sets of collisions [14,21,25], i,j,k,l € A={l,---,p}, and the 
summation is taken over all j,k,l, where A^j^^ arc nonnegative constants (related to transitional 
rates) satisfying A^'^ = ^^^^ = A^^, (u^ + u^- - - U;) = 0, A^-'^ = A'lj [14,25]. Fi is the 
(discrete) approximation of the collisional integral in the conventional continuous kinetic theory. 
The conditions defined for the discrete velocity above requires that elastic, binary collisions, 
such that momentum and energy are preserved. 

The collision operator is now simply obtained by joining A^j to the corresponding transition 
probability densities afj through A^j =S\ui-Uj \ afj, where, afj > 0, a^^j = 1, Ef=i = 1> 

V?, j = 1, • • • ,p; with S being the effective scattering or collisional cross-section [14,21,25]. If all 
n {p = 2n) outputs are assumed to be equally probable, then a^j=l/n for all k and I, otherwise 
a^j= 0. The term S'juj — Uj\dt is the volume spanned by the particle with Uj in the relative 
motion w.r.t. the molecule with Uj in the time interval dt. Therefore, S\ui —Uj\Nj is the number 
of j-particles involved by the collision in unit time. Collisions which satisfy the conservation 
and reversibility conditions which have been stated above are defined as admissible collisions 
[14,25]. 

Moreover, all the velocity directions after collisions are assumed to be equally probable. We 
note that, the summation of Ni {J2i^i) '■ the total discrete number density here is related to 
the macroscopic density : p{= nipY^iNi), where nip is the mass of the particle [21]. 
With the introducing of the Uehling-Uhlenbeck collision term [21] in equation (1) (Fj being 
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replaced or modified), 

= E WkNiil + jNi){l + jNj) - NiNj{l + jNk){l + jNi)], (2) 

for 7 < (normally, 7 = — 1) we obtain a gas of Fermi-particles; for 7 > (normally, 7 = 1) 
we obtain a gas of Bose-particles, and for 7 = we obtain equation (1) which is for a gas of 
Boltzmann-particles [14,24-25]. 

Considering binary collision only, from equation (2), the model of quantum discrete kinetic 
equation for Bose gases is then a system of 2n(= p) semilinear partial differential equations of 
the hyperbolic type (in two dimensional form): 

Wt^' + • ^iV, = ^ 5: iV,iV,-+„(l + 7iV,-+i)(l + 7AG-+„+i)- 

i=i 

2cSNiNi+n{l+lNi+i){l + ^Ni+n+i), i = l,---,2n, (3) 

where Ni = Ni^2n are unknown functions, and Uj =c(cos[6' + (i — l)7r/n], sin[0 + ~ 1)^/^]) 
are the particles velocities in the centre-of-mass coordinate system ; c is a reference velocity 
modulus and the same order of magnitude as the sound speed in the absence of scatters in [13] 
or [19], 9 is the orientation starting from the positive x— axis to the ui direction and could be 
thought of as a disorder induced by high pressure or external field (schematically shown in Fig. 
3). For example, there are admissible collisions (Ui,U3) < — > (U2,U4) as n=2 [25]. 
Since passage of the sound wave causes a small departure from an equilibrium state resulting in 
energy loss owing to internal friction and heat conduction, we linearize above equations around a 
uniform equilibrium state (particles' number density : Nq) by setting Ni{t,x) =No {1 + Pi{t,x)), 
where Pi is a small perturbation. The equilibrium state here is presumed to be the same as in 
Refs. [21-22]. After some similar manipulations as mentioned in Refs. [30], with B = jNq > 
(cf. Chu in [21], B defines the proportional contribution from the Bose gases; if 7 > 0, e.g., 
7 = 1), we then have 

where = (P^ -|- P^+„)/2, m = 1, • • • , n, since Di = for 1 = m (mod 2n). 
We start to look for the solutions in the form of plane wave -Dm= o,m exp i{kx — ut), (m = 
1, • • ■ , n), with uj=u{k) because it is related to the dispersion relations of quasi one-dimensional 
plane waves propagating in (monatomic) hard-sphere Bose gases. So we have 

(1 + ih(l + B)- 2\^cos\e + _ ^K^ + B) y „^ = q, m = 1, • • • , n, (5) 

where 

X = kc/{y/2u}), h{l + B) = hb = AcSNo{l + B)/u oc l/K^, 

h is the rarefaction parameter of the gas; Kn is the Knudsen number which is defined as the 

ratio of the mean free path of Bose gases to the wave length of the plane (sound) wave . 

Let am = C/{1 + ihb — 2A^cos^[^ + {m — l)7r/n]), where C is an arbitrary, unknown constant. 
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since we here only have interest in the eigenvalues of above relation. The eigenvalue problems 
for different 2 x n-velocity model reduces to 

3 Results and Discussions 

We can resolve the complex roots (A = \r+ i Aj) from the polynomial equation above and 
use the numerical way for direct plots. The roots are the values for the nondimensionalizcd 
dispersion (positive real part; a ratio of the sound or phase speed with respect to its continuum or 
hydrodynamical limit) and the attenuation or absorption (positive imaginary part), respectively. 
B could be related to the occupation number. We plot the main results into figures 4, 5, 6, and 
7, respectively. We firstly review the general characteristic dispersion relations for Bose gases 
before we interpret our present results. 

Curves in figures 4 and 5 (0 7^ 7r/4) resemble the conventional dispersion relations of ultrasound 
propagating in hard-sphere Boltzmann-statistic gases [30]. Present results show that as B or 
9 increases, the dispersion (A^ = krc/{-\/2uj)) will reach the hydrodynamical limit (as /i S> 1 
the wave speed is independent on the S or the s-wave scattering length, this result agrees 
qualitatively with Cowell et aZ.'s results [26] (for S > 0) or Bruun and Burnett's [28] (for 
B < 0)) earlier. That is to say, the phase speed of the plane wave in Bose gases (even for small 
but fixed h) increases w.r.t. to the continuum conditions {h 00) as the relevant parameter B 
increases or 9 increases (up to 7r/4). 

Meanwhile, as illustrated in Fig. 5, there always exist peaks or maximua in our attenuation 
results (related to the damping of the propagating wave). This agrees qualitatively with that 
of Jackson and Zaremba's result (for m=0 mode) [28]. Here, considering the Pauli-blocking 
effect, the maximum or peak absorption (or attenuation Aj = kic/ {\/2uj)) for all the rarefaction 
parameters h keeps the same for all B (say, B = 0.3 and B = 0.7). There are only shifts of the 
maximum absorption state (defined as h„iax) when B increases. It seems for the same mean free 
path (h oc the inverse of Kn) or mean collision frequency of the hard-sphere gases (i.e. the same 
h but h < hmax) there will be more absorption in larger 9 cases than those of = states when 
the plane wave propagates. 

To apply the acoustic analog, we should now introduce the brief idea [14]. In fact, studies of 
classical wave mechanical systems have some important advantages over quantum mechanical 
wave systems even there are similarities in-between [13]. In a mesoscopic system, where the 
sample size is smaller than the mean free path for elastic scattering, it is satisfactory for a 
one-electron model to solve the time-independent Schrodinger equation : 

-T^V^ip + V'(fU = Ei/j 
2m 

or (after dividing by — ^^/2m) 
+ [g^ - V{r)]^p = 0, 
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where q is an (energy) eigenvalue parameter, which for the quantum-mechanic system is \jlmE jfii^ . 
Meanwhile, the equation for classical (scalar) waves is 

or after applying a Fourier transform in time and contriving a system where c (the wave speed) 
varies with position r 

+ {(^ - y{f)\i) = 0, (7) 

here, the eigenvalue parameter g is w/cq, where a; is a natural (or an eigen-) frequency and cq is 
a reference wave speed. Comparing the time dependencies one sees the quantum and classical 
relation E = fiuj [14]. Thus, the energy E considered, e.g., in Refs. [19-20] for the resonant 
localization in quantum system corresponds to the frequency io considered here especially for 
the resonant states [14] by using the acoustical analogy. 

Back to our present figures w.r.t. /i, as now the parameter of disorder or free-orientation Q 
dominates. We can observe that, there is a continuous trend as Q increases toward 7r/4 ~ 0.7854. 
The dispersion (A^; a relative measure of the sound or phase speed) keeps increasing while the 
attenuation or absorption (Aj) keeps decreasing as Q increases from 0. At = 7r/4, there is no 
attenuation and dispersion, i.e., A^ = 1.0 and Aj = 0.0 [30]. This result also provides a good 
verification for the experimental side mentioned in [16-18] as there is no loss for this particular 
case (^ being a disorder parameter but fixed as 7r/4). We note that around the peak Aj state 
{hmin) as shown in Fig. 5, there exists a trend for the absence of diffusion (Aj starts decreasing 
rapidly) [18-20]. 

Based on the acoustic analog [14], from the definition of h or Kn, h = f collision / fsmmd-, where 
f sound is related to the classical frequency uj as mentioned above (cf. equation (7)) so that it is 
relevant to the energy E as defined for the localization, thus we can estimate the localization 
length from those figures which vary with h] f collision is the mean collision frequency of the 
particles. The localization length defined in [19] is proportional to the (hydrodynamic) mean 
free path / {I also depends on the internal frequency E as shown in [19]) and, comparing the 
definition of h here, is thus related to the inverse of h (or, say, the frequency) we used. We 
remind the readers that the temperature (T) could be linked to the mean free path or mean 
collision frequency under prescribed conditions [24] and the energy (thus to) could be related 
to the temperature (T) with the introdution of the Boltzmann constant kB- Based on these 
considerations, the relation for the localization length versus the frequency extracted from our 
results (especially in Fig. 5; the attenuation or absorption defined here is related to the inverse 
measure of, say, one wave length; the maximum absorption then corresponds to the minimum 
localization length in Fig. 5(a) of [19]) is qualitatively similar to that reported in [19]. This 
observation is now illustrated in figure 6 where we schematically define the localization length 
as the inverse of the wave absorption : l/Aj (per unit wave length). There also shows the 
exponential decay of the the localization length w.r.t. h or the inverse of the frequency (a 
corresponding measure to energy E in quantum- mechanic sense as already explained before). 
Thus we can also obtain similar results which resemble that reported in Fig. 5 (a) in [19]. 
We note that, as the rarefaction parameter {h = AcSNq/u) is fixed (cf. figure 6 for the schematic 
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localization length) , from its definition and the correspondence between the energy (E) and the 
frequency (u), the product of an effective scattering cross-section (S) (or the s-wave scattering 
length) and the number density of particles (Nq) must be proportional to £^ or a; once c and 
the localized state are made specific (the Pauli-blocking and disorder effects being excluded). 
Meanwhile, once or a; is fixed for the same situation mentioned above, the localized state will 
strongly depend on the S and Nq as well as B and 6. 

People might argue that a nonzero 6 would only make the system anisotropic, but not disordered. 
We should remind them that the derivation of present quantum kinetic approach was based on 
the binary encounter of a system of Bose particles. Once the mean free path and the centre of 
mass coordinate system were introduced (especially when the effective, admissible collision and 
the microreversibility which neglects the history and the correlations when particles traverse 
in phase space [14,21] were presumed) the randomness and disorder will occur although they 
are illustrated implicitly. The results presented here, in fact, also qualitatively resemble that 
reported in [29] where it was shown that when a periodic medium with a gap (in resulting 
spectra) is (slightly) randomized (like our disordered case : ^ / 0), possible localization occurs 
in a vicinity of the edges of the gap (like that of ^ = 7r/4 here; 6 = is implicit) [29]. As we 
only consider plane waves propagating in a hard-sphere gas, which is a kind of hard (Neumann) 
scatters [29], then it is interesting that our results for the dispersion relation resemble those of 
Neumann cases (especially Fig. 9) in [29]. 

To demonstrate the possible resonant states, we summarize our results in figure 7. Our results 
show that as 9 increases, the maximum absorption (A^) will decreases continuously except at 
9 = 7r/4 where there is a sudden jump (maximum)! This unusual absorption (or attenuation) 
peak is similar to those observations found in the Tx transition (temperature) for liquid helium 
(Bose liquids) (cf. Figs. 20, 21 in [31]). The interesting result is that, the absorption value 
obtained for B = 0.5 (Bose gases) is almost twice of that value for B = —0.5 (Fermi gases) at 
9 = 7r/4. Note that, once there arc Cooper pairings in Fermi gases, we can treat them (atomic 
pairs) as bosonic particles although the number density of them might be one half of the original. 
Thus, this resonant state might be relevant to the superfluid phase transition or Cooper pair 
formation tuned by the disorder. Meanwhile, once we consider the disorder is induced by the 
magnetic field, this sudden jump in Fig. 7 resembles that of resonance reported in [32]. 
As wc know, when a gas containing many identical particles is confined and cooled, the average 
momentum can be lowered so far that the typical de Broglie wavelength is larger than the average 
separation between the particles. In this case, the gas is said to be 'degenerate', meaning that 
the wave functions of neighboring particles overlap. Degenerate gases exhibit two dramatically 
different types of behavior, depending on whether the identical particles are bosons (such as 
photons) or fermions (such as electrons). 

To conclude in brief, our calculations here are useful to the understanding of waves propagating 
in microscopically random, disordered or granular media under strong confinement [2-3,5-7,9-12] 
and the dispersion management for atomic matter waves [16]. The direct relation of our results 
to the conventional one of static localization is not straightforward but could be understood 
qualitatively after the application of the acoustical-analog (i.e. the necessary transform from 
our complex A or real h to the conventional lo, E, V (potential barrier), and the characteristic 
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lengths (cf. [13-14,19]) : mean free path, wave length, etc.). Our results show that, as illustrated 
in figure 6, the resonant states strongly depend on the energy {E or u), the effective scattering 
cross-section (S), the number density (Nq) once the disorder and Pauli-blocking parameters are 
selected. We shall study other complicated problems in the future [2-3,33-35]. Acknowledgements. 
The author is partially supported by the China Post-Dr. Science Foundation (Grant No. 1999-17). 
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4 Appendix : Derivation of Eqn. (4) 

Prom Eqn. (3), after the linearization, we then have, (say, i=m) 

— Pm + ■ g^Pm + 2cSNq[{P^ + P^+n) + ^NoiPm + Pm+n + Pj^) + ■ ■ ■] = 

^ ^[(Pfe + Pk+n + jNoiPk + Pk+n + Py) + -- ■], (8) 

here, m = 1, • • • , 2n, Pj^ = for n=2 because of the restriction for the total perturbations in 
an equilibrium state and the remaining terms in both sides are higher order terms related to 
(7A^o)^- The linearized version of above equation (for n=2) is 

o o O C AT 

-Pm + Vm ■ g^Pm + 2cSNoiPm + Pm+n)(l + T^^o) = E ^^{1 + jNo). (9) 

In these equations after replacing the index m with m + n and using the identities Pm+2n = Prm 
then we have 

Q-^Pm+n - " Q^Pm+n + 2cSNo{Pm + Pm+n)(l + TA^o) = E ^^(1 + TA^o)- (10) 

Combining above two equations, firstly adding then subtracting, with A^, = {Pm + Pm+n)/2 
and Bm = {Pm - Pm+n)/2, we can have 

^Am - c cos[e + + 4cSNoAm{l + 7iVo) = E + 7iVo), (11) 

OT n OX n f-^^ 

k=l 

—Bm + ccos[e+ ^'^~^ '^ ]—Am = 0, m = 1, • • • , 2n. (12) 
From Pm+2n = Pm, and with Am = (Pm + Pm+n)/2 and Bm = {Pm - Pm+n)/2, we can have 

After eliminating Bm, with B = jNq, we then have 

2 2r« {m - l)7r , . r..r r..d,^ AcSNo{l + B) d ^ 

where Dm = {Pm + Pm+n)/2, m = 1, ■ ■ ■ ,n, since Di = Dm for 1 = m (mod 2n). 



^4 



^4 



Fig. 3 Schematic plot for the regular scattering and the disorder-influenced scattering. 
Plane waves propagate along the X-dircction. Binary encounters of Ui and C/3 and their 
departures after head-on collisions {U2 and U4). Number densities Ni are associated to Ui 
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Fig. 4 Disorder or orientational (0) effects on the dispersion (A^). 
h = 4cSNq/lj, S is the effective scattering cross section, Nq is the 
number density, B = jNo is the Pauli-blocking parameter. 
This result agrees with Cowell et a/.'s result [26] for the hydrodynamic 
regime {h ~ 0(10) here). Wave speed is independent on S for larger h. 
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Fig. 5 Disorder or orientational (0) effects on tfie attenuation (Aj). 
Tfiis result agrees witfi Jackson and Zaremba's result (m=0 mode) [28]. 
Tfiere always exists a peak or maximum A^. 
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Fig. 6 Disorder {9) effects on tfie localization lengtfi (l/Aj). h = AcSNq/uj. 
Note tliat tlie energy E corresponds to tiuj [13-14,19]. This figure is a scfiematic type, 
(of. the presentation : Fig. 5 (a) for that used in [19].) All units are dimensionless. 
As B increases, the minimum (for the localization length) state hmin decreases, i.e., 
the temperature becomes much more lower (as the mean free path is rather large). 
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Fig. 7 Variations of (max.) Aj w.r.t. the disorder or free orientation : for 

B= ±0.5. The sudden jump at ^ = 7r/4 imphes the resonant transition in Bose fluids 

may also occur for Fermi fluids (in bound pairs). Aj for B = 0.5 (Bose gases) 

is almost twice of that value for B = —0.5 (Fermi gases) at ^ = 7r/4. 



